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Abstract 

We study the impact of recent 5-factories measurements and upper limit of ra- 
diative and semileptonic rare -B-decays. We present model independent constraints 
on the relevant Wilson coefficients and show the impact on the parameter space of 
some concrete realizations of the Minimal Supersymmetric Standard Model. 



1 Experimental inputs 



In this talk we use the most recent experimental results on inclusive and exclusive b — > 57 
and b — > s£ + £~ decays to update the analysis presented in Ref. ||T|. The experimental 
results that we use in the analysis are 

B{B^X sl ) = (3.4018$ ) x ID" 4 (1) 

B(B^K£ + r) = (0.63^) x 10- 6 1J, (2) 

B{B -> K*n + /T) < 3.0 x 10~ 6 at 90% C.L. §J, (3) 

B(B -> K*e + e~) = (1.68±g;|| ± 0.28) x 10~ 6 §, (4) 

B(B^X sf i + f,-) = (7.9±2.1^;°)xl0- 6 @, (5) 

B(B -> X s e + e~) = (5.0 ± 2.3±\$) x 10~ 6 @, (6) 

B(B -> x s e + r) = (6.1 ± lAtU) x 10~ 6 §. (7) 



The B — > X s 7 branching ratio is given with a cut on the photon energy {E 1 > m&/20) 
and the value presented is the weighted average of the four available measurements. The 
branching ratios for the semileptonic modes (0-0) refer to the non-resonant branching 
ratios integrated over the dilepton invariant mass spectrum. On the theoretical side, 
this amounts to consider only the perturbative part of the amplitude and to leave out 
all the resonant contributions (that are usually shaped using Breit-Wigner ansatze). On 
the experimental one, judicious cuts are used to remove the dominant resonant contribu- 
tions arising from intermediate cc resonances (J /if), ip', ...); SM-based theoretical distribu- 
tions |J are then used to correct the data for the experimental acceptance. Finally, let us 
stress that the B — > X s e + e~ branching ratio is given with a cut on the di-lepton invariant 
mass, M ee = y^i > 0.2 GeV, in order to remove virtual photon contributions and the 
7t° — > ec-f photon conversion background. The B —> X s e + e~ rate increases steeply for 
s — » (because of the almost real photon pole) and is extremely sensitive to the Wilson 
coefficient of the magnetic moment operator (that controls the b — > 57 transition). The 
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use of the branching ratio without extrapolation to the full s spectrum, allows to reduce 
the uncertainties in the comparison with the Standard Model (SM) prediction and to 
properly analyze models with an enhanced magnetic moment Wilson coefficient. 



2 Theoretical framework and SM predictions 

The effective Hamiltonian in the SM inducing the b — > s£ + £~ and b — > 57 transitions is 
(up to negligible contributions proportional to V* s V ub ) 
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where Gf is the Fermi constant and are the CKM matrix elements. Oj(/i) are 
dimension-six operators at the scale /z and Cj(/i) are the corresponding Wilson coeffi- 
cients. The most relevant operators are (the complete list can be found, for instance, 
in Ref. M) 
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where the subscripts L and R refer to left- and right- handed components of the fermion 
fields. 

The differential decay width for the inclusive decay B — > X s £ + £~ is given by the parton 
level result supplement by calculable power corrections. In the NNLO approximation, the 
non-resonant decay width can be written as 
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where s = s/ml and the functions Gi(s) (i = 1,2,3) and G c (s) encode respectively the 
l/m-l and corrections. Cf^ are effective Wilson coefficients (whose explicit form 

is given in Ref. 0) that are functions of the dilepton mass squared and incorporate 

part of the operator matrix elements. In particular, Cf^ contains contributions due to 
perturbative cc rescattering and develops an imaginary part for s > Am? c . Performing the 



s integration and constraining s ee > (0.2 GeV) 2 , the decay widths in electrons and muons 
are essentially equal and are given by the following numerical formula: 



B(B -> XJtiT) = [4.534 + 8.665 |C* ot | 2 + .119 (|C 9 NP | 2 + |Cf p | 2 ) + .996 ReC* ot C, 
+4.130 ReC* ot + 0.171 ImCf 1 + 1.068 ReC 9 NP + .064 ImC 9 NP -1.011 ReC^l x 10 
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where Cg P and Cf p are the new physics contributions to C 9 (/iy^) an d C w (fiw) evaluated 
at fi w ~ m w and C^ ot is the sum of the SM (Cj M (/j,b)) and new physics (Cf p (fib)) 
contributions evaluated at fib — 2.5 GeV. A detailed discussion of all the assumptions 
that enter this formula can be found in Ref. [|lj]. In the SM we find 

B(B -> XJ+tr) = (4.15 ± 0.27 ± 0.21 ± 0.62) x 10~ 6 = (4.15 ± 0.70) x 10~ 6 (14) 

where the errors correspond to variations of fib, nit,poie an d m c /mb. Comparing this 
estimate with the experimental measurements ([§)— (|7|) we see that there is agreement 
with the SM at the 1 a level. Note that Previous experimental data on the di-electron 
channel were extrapolated using SM-based distributions and that the branching ratio for 
B — > X s e + e~ integrated over the whole dilepton invariant mass spectrum is |l] (6.89 ± 
1.01) x 10 -6 (i.e. the s ee < (0.2 GeV) 2 region enhances ([14]) by 66%). This clearly shows 
to what extent the choice to give branching ratios not extrapolated allows for a cleaner 
identification of new physics effects. 

For what concerns the exclusive decays B — > K(*'£ + £~, we implement the NNLO 
corrections calculated by Bobeth et al. in Ref. 0] and by Asatrian et al. in Ref. [jlTJ for 



the short-distance contribution. Then, we use the form factors calculated with the help 



of the QCD sum rules in Ref. |L2|]. For lack of space we can not describe some subtleties 
related to the treatment of the so-called hard spectator interactions and to the value of 
the magnetic moment form factor at s = (see Ref. |1[ for a complete discussion). The 
SM NNLO predictions that we obtain are 

B(B -> KtT) = (0.35 ± 0.12) x 10~ 6 (5B KU = ±34%) , (15) 
B(B -> K*e + e-) = (1.58 ± 0.49) x 10~ 6 (5B K * ee = ±31%) , (16) 
B(B -> K*fi + fi-) = (1.19 ± 0.39) x 10~ 6 (8B K *^ = ±33%) . (17) 

From the comparison with Eqs. ©-(D) we see that in the B — > Kt + i~ channel there is 
a 1.6 a discrepancy with the SM expectation while in the K* channels there is a perfect 
agreement. 



3 Model independent analysis 

The first step consists in extracting the bounds that the measurement ([!]) implies for 
C7 0t (2.5 GeV). The main difficulty arises from the treatment of the m c dependence of the 
B — > Xg'y branching ratio. In Ref. [0, it was noted that, in this decay, the charm quark 



mass enters the matrix elements at the two-loop level only and that it would be more 
appropriate to use the running charm mass evaluated at the fib — 0{m b ) scale, leading 




Figure 1: 90% C.L. bounds in the [R 7 (fi), R 8 (fi)} plane following from the world average 
B — ► X s j branching ratio for \i = mw (left-hand plot) and fi = 2.5 GeV (right-hand plot). 
Theoretical uncertainties are taken into account. The solid and dashed lines correspond 
to the m c = m^poie and m c = m^ IS (fj J b) cases respectively. The scatter points correspond 
to the expectation in MFV models. 



to m c /mfc ~ 0.22 ± 0.04, compared to m CjPO ; e / 'mi, ~ 0.29 ± 0.02. This is a reasonable 
choice since the charm quark enters only as virtual particle running inside loops; formally, 
on the other hand, it is also clear that the difference between the results obtained by 
interpreting m c as the pole mass or the running mass is formally a NNLO effect. In what 
concerns b — > s£ + £~ , the situation is somewhat different, as the charm quark mass enters 
in this case also in some one-loop matrix elements. In these one-loop contributions, m, 



c 



has the meaning of the pole mass when using the expressions derived in Ref. ||11|| . Since 
the bounds on the Cj do not depend dramatically on m c , we just derive them using both 
values of the charm mass and taking the union of the allowed ranges. We present the 
results of this analysis in Figs. |TJa and [I|b, where we show the allowed regions in the R7 
and R 8 plane obtained using the 90% C.L. B -> X sl bound (here R 7 , 8 = C^/Cff). 
We take \R 8 (fj,w)\ < 10 in order to satisfy the constraints from the decays b — » sg and 
B — > X(f [|14]].The regions in Fig. |l|b translate in the following allowed constraints: 

f m c /m b = 0.29 : C* ot (2.5 GeV) G [-0.37, -0.18] & [0.21, 0.40] , , , 

\ m c /m h = 0.22 : C\° l (2.5 GeV) G [-0.35, -0.17] & [0.25, 0.43] . 1 ' 

In the subsequent numerical analysis we impose the union of the above allowed ranges 

- 0.37 < C7 Ot ' <0 (2.5 GeV) < -0.17 & 0.21 < C7° t,>0 (2.5 GeV) < 0.43 (19) 

calling them C^-positive and C^-negative solutions. 

We present the results of the model independent analysis of b — > s£ + £~ decays in Fig. [2| 
Within each plot we vary C 7 ot inside the allowed ranges ( |19D and plot the 90% C.L. 
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Figure 2: NNLO Case. Constraints on new physics contributions to the Wilson co- 
efficients Cg and Cio implied by b — > s£ + £~ decays. The plots correspond to the 
C7 0t (2.5 GeV) < and Cj ot (2.5 GeV) > case, respectively. The points are obtained 
by means of a scanning over the EMFV parameter space and requiring the experimental 
bound from B — > X s 7 to be satisfied. 



constraints implied by Eqs. (0)-(0) in the [Cg P (fiw) , Cf^} plane. The SM correspond 
to the point (0,0). In each plot the inner and outer contours are determined by the 
measurements of the decays B — > K£ + £~ and B — > X s £ + £~ respectively. 



4 Analysis in supersymmetry 

In this section we analyze the impact that the measurements have on three vari- 

ants of the minimal supersymmetric standard model (MSSM), namely minimal flavour 
violation (MFV), gluino mediated contributions and extended minimal flavour violation 
(EMFV). 

MFV. As already known from the existing literature (see for instance Ref. ||15||), 
minimal flavour violating contributions are generally too small to produce sizable effects 
on the Wilson coefficients Cg and C\q. Indeed, scanning over the MFV parameter space 
and imposing the lower bounds on the sparticle masses we obtain 

r tot < n . f C 9 MFV (^ W ) G [-0.2, 0.4] tot f Cl IFV ^ w ) G [-0.2, 0.3] , , . 

° 7 <u - \C^ FV e [-1.0,0.7] ' ° 7 >u ' \C*< FV e [-0.8,0.5] ' 1 Uj 

From the comparison of the size of these contributions with the allowed regions depicted 
in Fig. ^] we see that the current experimental results on b — > s£ + £~ decays are not precise 
enough to constraint the MFV parameter space. The situation is completely different for 
what concerns b — > sj. The scatter plot presented in Fig. [l] is obtained varying the MFV 
SUSY parameters and shows the strong correlation between the values of the Wilson 




Figure 3: Left plot: Dependence of Rf (fib) = (jib) I 'Cf M (fib) on the mass of the 
charged Higgs. Right plot: Dependence of R* (fib) = C*(fib) /C^ M (fib) on the mass of 
the lightest stop in MFV models. The chargino contribution is essentially proportional to 
sin 9f tan (3 for not too small sin 6^-. For the curve 2 we show the variation due to several 
choices of 9^ and tan/?. 



coefficients CV and C%. In fact, the SUSY contributions to the magnetic and chromo- 
magnetic coefficients differ only because of colour factors and loop-functions. In Fig. || 
we present, finally, the dependence of the charged Higgs and chargino contributions to 
C-j on the relevant mass parameters (that are the charged Higgs mass for the former and 
the lightest chargino and stop masses for the latter). From these figures is clear that 
the knowledge of the sign of C^ ot will strongly constrain the MFV parameter space, by 
putting upper or lower limits on the chargino and stop masses. 

Gluino contributions. Gluino contributions to Cg and Cio are governed by mass 
insertions in the down squark mass matrix. From the analysis presented in Ref. |15| we 
see that the dominant diagrams involve the parameter (S^ll and that large deviations 
from the SM are unlikely. 

Extended— MFV models. A basically different scenario arises if chargino-mediated 
penguin and box diagrams are considered. As can be inferred by Table 4 in Ref. JHJ, 
the presence of a light t 2 generally gives rise to large contributions to Cg and especially 
to Ciq. EMFV models are based on the heavy squarks and gluino assumption. In this 
framework, the charged Higgs and the lightest chargino and stop masses are required to 
be heavier than 100 GeV in order to satisfy the lower bounds from direct searches. The 
rest of the SUSY spectrum is assumed to be almost degenerate and heavier than 1 TeV. 
The lightest stop is almost right-handed and the stop mixing angle (which parameterizes 
the amount of the left-handed stop II present in the lighter mass eigenstate) turns out 
to be of order 0(mw/Mg) ~ 10%. The assumption of a heavy (> 1 TeV) gluino totally 
suppresses any possible gluino-mediated SUSY contribution to low energy observables. 
Note that even in the presence of a light gluino (i.e. M- g ~ 0(300 GeV)) these penguin 
diagrams remain suppressed due to the heavy down squarks present in the loop. In the 
MIA approach, a diagram can contribute sizeably only if the inserted mass insertions 
involve the light stop. All the other diagrams require necessarily a loop with at least 
two heavy (> 1 TeV) squarks and are therefore automatically suppressed. This leaves 



us with only two unsuppressed flavour changing sources other than the CKM matrix, 
namely the mixings u L — i 2 (denoted by 6q $ ) and cl — i 2 (denoted by S 5 j 2 ). We note 
that 8q 2 an d ^c L t 2 are mass insertions extracted from the up-squarks mass matrix after 
the diagonalization of the stop system and are therefore linear combinations of (#13)^, 
(£l3)ll and of (623)^1 (^23)^5 respectively In Fig. |2] we present the results of an high 
statistic scanning over the EMFV parameter space requiring each point to survive the 
constraints coming from the sparticle masses lower bounds and b — *> S7. Note that these 
SUSY models can account only for a small part of the region allowed by the model 
independent analysis of current data. In the numerical analysis reported here, we have 
used the integrated branching ratios alone to put constraints on the effective coefficients. 
This procedure allows multiple solutions, which can be disentangled from each other only 
with the help of both the dilepton mass spectra and the forward-backward asymmetries. 
Only such measurements would allow us to determine the exact values and signs of the 
Wilson coefficients C 7 , C 9 and C w . 
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